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Abstract 

We study the formation of low-mass black hole X-ray binaries with main sequence 
companions that have formed through case C mass transfer (mass transfer following 
the helium core burning phase of the black hole progenitor). We identify these 
objects with the observed soft X-ray transients. Although this scenario requires a 
set of fine tuned conditions, we are able to produce a current Galactic population 
of ~2000 objects, in agreement with estimates based on observations. 

The narrow interval in initial separations leading to case C mass transfer, com- 
bined with the allowed narrow range of separations after the common envelope evo- 
lution, constrains the common envelope efficiency in this scenario: Xuce ~ 0.2 — 0.5. 

Key words: black hole physics — stars: binaries: close — accretion 



E-mail: popenoe@nuclear.physics.sunysb.edu 
^ E-mail: chlee@kias.re.kr 
^ E-mail: tauris@nordita.dk 



Preprint submitted to Elsevier Preprint 



5 February 2008 



1 Introduction 



In the literature there have been severe difficulties in evolving a sufficient num- 
ber of black hole soft X-ray transient sources (SXTs) Ergma & van den Hcuvel 
(1998), Portegies Zwart, Verbunt, & Ergma (1997). For example, Portegies 
Zwart et al. obtain a birth rate of 9.6 x 10~^ yr~^ for the SXTs, to be com- 
pared with one of 2.2 x 10^^ yr~^ for the binaries with a neutron star, whereas 
they remark that observations indicate equal formation rates for these two 
types of binaries. Their large discrepancy arises from their limit ZAMS mass 
40Mq for black hole formation. 

In fact, however, for the evolution of the SXTs, a ZAMS mass more like 
~ 20Mq should be chosen as limit for high-mass black hole formation. This 
is the limit for single stars (Brown ct al. 2001). If the massive black hole 
progenitor in the SXTs can complete He core burning before it is removed in 
common envelope evolution (Case C mass transfer) then it will evolve like a 
single star since the remaining ~ 10^ yrs of its hfetime is too short for He wind 
losses to effect its evolution. Indeed, Ergma & van den Heuvel (1998) realized 
that '^Mbh should not be larger than 20 to 25M0" in order to have sufficient 
SXTs, but they did not realize that the more massive stars in a binary evolve 
in a different way from a single massive star. 

If the H envelope of the more massive star in a binary is removed by early 
Roche Lobe overflow. Case A or Case B, before He core burning is completed, 
the resulting "naked" He core will blow away, leaving an Fc core too low in 
mass to evolve into a high- mass black hole (Brown et al. 2001). 

The common envelope evolution of Portegies Zwart et al. (1997) for a 2OM0 
turns out to be useful for us, because this will be the same, regardless of 
whether the resulting He star goes into a neutron star as they believed, or into 
a high-mass black hole. It turns out that this evolution can only be carried 
out in Case C, so the He star is clothed during He core burning, appropriate 
for our high mass black hole evolution. 

Given the Schaller et al. (1992) evolution, we find that Case C evolution does 
not work for stars more massive than ~ 22Mq, in agreement with Portegies 
Zwart et al.(1997). However, we need the more massive black hole progenitors 
to evolve the transient sources with subgiant companions, as we shall discuss. 

In this note we consider the evolution of those low-mass black hole X-ray bi- 
naries with main sequence companions which are formed through case C mass 

transfer followed by common envelope and spiral-in evolution. We estimate 
their expected number currently present in our Galaxy. We therefore discuss 
the following. 
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(1) The initial separation interval of binaries that go through case C mass 
transfer and the ZAMS mass interval of stars being able to produce high- 
mass black holes (Sect. 2). 

(2) The initial masses of the black hole companions (the donors in the present- 
day X-ray binaries) and the allowed orbital periods after the common 
envelope required for the systems to become observable X-ray binaries 
within the age of our Galaxy before the donors leave the main-sequence 
(Sect. 3). 

(3) The resulting common envelope efficiency that can be estimated combin- 
ing points 1 and 2 (Sect. 4). 

(4) The evolution of the mass-transfer rate in the X-ray transient phase and 
the lifetimes of the systems (Sect. 5). 

(5) The total number of systems expected in our Galaxy at present (Sect. 6). 

Following the discussion of the above 5 points, wc briefly discuss the black 
hole systems with (sub)giant companions Nova Scorpii (GRO J 1655— 40) with 
a 2AMq F-star companion and a 6.3M0 black hole and V4641 Sgr (XTE 
J1819— 254) with a 6.5M0 B-star companion and a Q.GM© black hole. At 
least in the former case there has been substantial mass loss in the black 
hole formation, SM©. Thus, a He core of at least IIM© is needed for the 
black hole progenitor, or a ZAMS 35Mq star. From this we conclude that 
the Schaller ct al. (1992) wind losses are too large and that they must be 
substantially decreased in order to evolve these two binaries. 



2 Case C mass transfer: Limits on Initial Sepeirations and Black 
Hole Progenitor Masses 

In order to have case C mass transfer the radius of the star has to expand after 
core helium burning has stopped. In the calculation of Schaller et al.(1992), 
stars with ZAMS masses 2OM0 and 25M0 satisfy the requirement. However, 
during the core helium burning stage, the orbit expands due to the mass loss: 

a M-AM + Md ^ ' 

where a is the orbital separations, M is the mass of the black hole progenitor, 
AM is the mass lost from the black hole progenitor during the core helium 
burning, and is the donor star mass. In order to initiate mass transfer after 
core helium burning the star has to expand sufficiently that this widening of 
the orbit is compensated for. 

In the upper curve of Fig. 1 we draw the orbital separation of the binary 
required to initiate Roche lobe overflow calculated at each stage from the 
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Fig. 1. Radius of black hole progenitors (R) and the initial orbital separations (oj) 
of the progenitors of X-ray transient binaries with a IMq companion. Upper dot 
dashed curves correspond to the orbital separations required to initiate Roche lobe 
overflow at the radius taken for each stage from Schaller et al. (1992). The solid 
curves correspond to the required initial separations after corrections of the orbit 
widening due to the wind mass loss. The lower dotted curves correspond to the 
radius of the black hole progenitors taken from Schaller et al. (1992). 

radii given by Schaller et al. (1992). For the masses considered the radius of 
giant star is at ~ 2/3 of the radial separation between the two stars. The 
solid curve gives the corresponding initial separation between the massive 
star and the low- mass main sequence companion after strong wind mass loss 
is switched on. Here we assume that the wind mass loss does not affect the 
stellar radius.[^ We see that in the case of a 20Mq star the mass transfer 
can start following He core burning for binaries with an initial separation in 
the interval of 137Oi?0 < < 16O5i?0. In the case of a 25Mq star, however, 
mass transfer is only possible up to just after the beginning of core helium 
burning (with initial separations around 1170 Rq). Binaries with larger initial 
separations will become too wide during core helium burning so even though 
the star expands slightly after core helium burning, it will not be able to fills 



^ One of us (T.M.T.) has carried out stellar evolution calculations, using Eggleton's 
code, with/without the wind mass losses, and find that the radius behaves almost 
the same. The minor difference will give rise to adjusting our numbers slightly, but 
the outline will remain the same. 
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its Roche lobe (see also Fig. 2 of Portegies Zwart, Verbunt, & Ergma 1997). 

Brown, Weingartner & Wijers (1996) found that stars with ZAMS masses 
above 19M0, evolved by Woosley & Weaver (1995), evolve into high-mass 
(~ 7Mq) black holes. As outhned in Brown et al. (2001), this requires the 
special rate of 170 keV barns at energy E = 300 kcV for the ^^C(a,7)^^0 
reaction used by Woosley & Weaver. Recent experiments including both El 
and E2 components, obtain Sf^^ = (165 ± 50) keV barns (Kunz et al. 2001). 
Brown et al. (2001) discuss quantitatively how changes in this rate change the 
ZAMS mass at which the calculated Fe core mass increases rapidly with mass. 
The fact that the progenitor mass of SN1987A lay in the narrow interval of 
18 — 2OM0 appropriate for forming black holes with the Woosley and Weaver 
170 keV barns gives further support to this value. A ZAMS ISM© star, which is 
usually taken as progenitor for SN 1987A, will form a low-mass (~ 1.5— l.SM©) 
black hole or a neutron star. 

Hence, the mass interval of stars that can go through case C mass transfer and 
collapse into a high-mass black hole is hmited between 19 and 20-25 Mq. We 
shall (arbitrarily) take 22 Mq as an upper limit. Stellar evolution calculations 
with lower wind mass losses might give a higher mass limit. 

Using the above derived limits on the mass of the black hole progenitors and 
the initial separation, the rate of formation of black hole transient sources with 
main sequence companions we can evolve (see Section 6) is nearly an order of 

magnitude less than in Brown, Lcc & Bethe (1999) where wind mass losses 
were ignored. We realize that our evolution requires a fine-tuned parameter 
space, but we shall argue that it all hangs together. The exact numbers may 
alter with changes in stellar evolution models. 



3 Limits on Donor Meisses and Orbital Periods 

During the common envelope evolution the low-mass main sequence star of 
mass ~ IM0 will spiral into the envelope of the massive giant. The orbit will 
shrink dramatically and the outcome is a close binary consisting of the core 
of the massive giant and the low-mass star (Paczyhski 1976). Shortly after 
the spiral-in and ejection of the envelope of the giant, the remaining core will 
collapse into a black hole. Whether or not this collapse is associated with 
mass loss is still an open question. The orbit will then shrink further due to 
the loss of orbital angular momentum via magnetic braking and gravitational 
wave radiation (e.g. Verbunt 1990) until the low-mass (donor) star begins to 
transfer mass to the black hole - forming an X-ray binary. Since the wind 
mass-loss rate of the low-mass donor star is very low, and its accretion during 
the short common envelope phase is negligible, it will have a mass at the onset 
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of the X-ray phase which is about equal to its initial mass. 



The time evolution of the orbital separation is described as 



a J orb J orb V ^V-lBH^v^d / 

Here we assume that the mass lost from the companion star is all accreted 
onto the black hole. The orbital angular momentum loss by gravitational wave 
radiation (gw) and magnetic braking (mb) is given by: 
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^ -0.5 X + ^'^^' s- (3) 

where is the radius of the donor star, hmb is an efficiency parameter given 
below, and 7 Ri 0.1 MdR^ is the moment of inertia of the donor star. Here we 
have assumed a magnetic braking law based on observations of slowly rotating 
single stars by Skumanich (1972). (However, this law may be too strong, i.e. 
overestimating the dependence on the angular velocity, see e.g. Stepien 1995). 
We further assumed the the magnetic braking efficiency suggested by Kalogera 
& Webbink (1998), 



bmbiMd) 



0, Md < O.37M0 

1, 0.37Mq <Md< 1.O3M0 . (4) 
exp(-4.15(Md - 1.03)), Md > IMMq 



Here the effect of magnetic braking is assumed to be strongest for stars with 
a mass ^ I.OMq. Donor stars with mass > 1.5M0 do not have a convective 
envelope on the main sequence and the magnetic braking effects are therefore 
negligible, so the binary evolution is dominated by gravitational radiation until 
the mass transfer starts. 

In Fig. 2, wc show the possible ranges of orbital periods and donor masses 
immediately after the formation of the black hole, such that an X-ray binary 
with a main-sequence donor is formed within the Hubble time (Kalogera 1999). 
If no mass is lost in the formation of the black hole, these periods are equal 
to the periods immediately after the common envelope as a result of the short 
time interval before the collapse of the core. 

The upper limit is set by two conditions: i) the binary should start mass 
transfer within the Hubble time (line denoted by A), ii) the mass transfer 
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Fig. 2. Limits on the orbital period after the formation of the black hole as func- 
tion of the mass of the low-mass companion. The upper and lower boundaries (thick 
solid lines) of mass and period are set by the conditions (see Section 3): A) Mass 
transfer begins within the Hubble time, B) Mass transfer begins before the com- 
panion evolves off the main sequence, C) Mass transfer from the companion begins 
right after finishing common envelope evolution. The expected ranges of periods 
and companion masses, after the common envelope phase, are plotted for various 
given efficiencies (see Section 4). The width of each band for an assumed Aace 
is determined by the initial band of possible Roche Lobe overflows with Case C 
mass transfer. We took the ZAMS mass of the black hole progenitor to be 20Mq 
which corresponds to Mp ^ 16. SM© in the beginning of case C mass transfer, and 
Mue = Mbh = 7Mq. 

should start before the donor evolves off the main-sequence (line denoted by 
B) which also determines the lower limit on Aace = 0.2 (see Section 4). The 
kink in line B around ^ IM© comes from the fact that the effect of magnetic 
braking, eq.(4), is strongest for stars with a mass ~ IMq. For this plot, we 
took the simple interpolation of Schaller's results for the evolution of a IMq 
star with two extrapolations for the radius R and the lifetimes of the main- 
sequence star tend ms; 



R{Md,t)=RM^{t) X M° 



0.88 












(5) 



where tju© is the time at the stage 13 of the IMq model of Schaller et al. 



1992 (sec Fig. 1). The applied formulae may not be very accurate, but the 
qualitative behaviour after more realistic numerical calculations will remain 
the same. The lower boundary of the period immediately after the common 
envelope (line denoted by C) is set by the condition that the main sequence 
star should not overfill its Roche-lobe at the end of the common envelope. 



4 Common Envelope Evolution: Binding Energy and Efficiency Pa- 
rameter 

Since our initial separation Oi is sharply defined by the condition of case C 
mass transfer, and our final separations are constrained by the distance that 
can be traversed during magnetic braking and gravitational wave radiation, 
we can determine the range of allowed common envelope efficiencies in order 
to form low-mass black hole X-ray transients. 

During the common envelope phase the energy needed to expel the hydrogen 
envelope of the black hole progenitor is tapped from the drop in binary orbital 
potential energy with efficiency ctce (also denoted r) in the literature): 

GMpMe _ GMpM^ _ fGMffeMd GMpM^ 
XR Xridi '^^ \ 2a f 2ai 

where Mp is the total mass of the BH progenitor star just before the common 
envelope forms, Mg is the mass of its hydrogen envelope, Mhc is the mass of 
its core, a\ and af is the initial and final separation, before and after the 
common envelope, respectively, tl is the dimensionless Roche-lobe radius. 
Given the parameters of the system at the start of the common envelope, 
the final separation is determined by A, describing the structure of the giant 
and dee, the efficiency of the energy conversion. In our case the final separation 
is limited as shown in Fig. 2 so the product Attce can be constrained. In the 
literature A = 0.5 had often been used, and a high efficiency ace > 1 was 
often required in order to explain the observations. However, recent detailed 
stellar evolution calculations by Dewi & Tauris (2000) show that A can be 
substantially larger. These high values of A solve the problem of unrealistically 
high efficiencies ace- For a 2OM0 star, corresponding to I6.2M0 at the tip of 
the AGB, Tauris & Dewi (2001) find A- values in a large interval: 0.1 < A < 3 
depending on the exact location of the core mass boundary and amount of 
internal thermodynamic energy included. 

In Fig. 2 we plotted the possible ranges of periods after the common envelope. 
We took the ZAMS mass of the black hole progenitor to be 2OM0, which 
corresponds to Mp ~ 16.8Mq in the beginning of case C mass transfer, and 
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Fig. 3. Time evolution of the donor star mass and mass loss rate for M^^i = 1.25Mq 
with a 7Mq black hole and an initial period of 20 hours. The dot-dashed hnes on the 
right panel correspond to the critical mass loss rate for the steady X-ray sources, 
eq. (8). The two tracks correspond to the different time scales of stripped donor 
stars (see the discussion of Sec. 5). 

assumed Mhc = Mbh = 7Mq. For each value of A«ce the two lines are for 
the limiting initial separations Oi (see Sect. 2). Donor stars with masses above 
I.SMq can only be formed if the common envelope efficiency Xace is around 
0.2. 



5 Life Time of X-ray Transient Sources 



Brown, Lee & Bethe (1999), assuming a mass-loss rate of 10~^Mq yr~^, ob- 
tained a lifetime of the X-ray transients of 10^ yrs . However, as Ergma & 
Fedorova (1998) discussed, the mass-loss rate changes as a function of time. 

Once the Roche lobe overflow starts, at any given time, the orbital separation 
is self determined in order for the companion star to fill its Roche lobe. 



a^Rd/rL{Md,MBH) (7) 

Hence, the mass loss rate is automatically determined by the feedback effects 
between the orbital widening due to the mass loss and the orbital contraction 
due to the gravitational wave radiation and the magnetic braking, eq. (2). In 
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the numerical simulation of mass transfer due to Roche lobe overflow, for a 
given 5Md, we can get the 5t, equivalently M^, by requiring the donor star fill 
its Roche lobe. 

In Fig. 3 we plot the evolution of the donor star mass and the mass loss rate 
as a function of time for M^ j = 1.25M0 with a 7Mq black hole. There are 
two tracks after Roche lobe overflow. 

(1) Case TO: we assumed that the stripped star after Roche lobe overflow 
continues to follow exactly the same evolution time scale as the original 
star. 

(2) Case TS: wc assumed that the stripped donor star loses all its initial 
information. When the star become smaller in mass, it follows the time 
scale of the main sequence star of that reduced mass. 

We beheve the realistic situation to he between case TO and TS. The radius of 
the donor star, which essentially determines the orbital separation, is smaller 
in case TS because the donor star is less evolved than that in TO case. Even 
though, for any given donor star mass M^, the lifetime of the donor star itself 
is longer in TS case, the smaller radius of the donor star gives the shorter 
orbital separation with larger mass loss rates, finally decreasing the lifetime 
of X-ray transients in case TS. 

In Fig. 3, the time scales of X-ray transients are ~ 2 x 10^ yrs for a 1.25M0 
donor star. The mass of the donor star (initially a 1.25Mq star) drops quickly 
in the beginning with a high mass loss rate > IO^^Mq yr"-*^, which slows down 
afterwards. Any donor stars > O.SM© show the similar behaviour. Hence, 
independently of initial donor masses, most of the time X-ray transients will 
appear as low-mass X-ray transients (i.e. with a donor star mass < Q.'oMq) 
having a mass-loss rate of ~ 1O~^°M0 yr~^, which is consistent with the 
empirical rate of 1.3 x IO^^^Mq yr~^ found by Van Paradijs (1996). The dot- 
dashed lines on the right panel of Fig. 3 correspond to the critical mass loss 
rate for the steady X-ray sources (King et al. 1997), 



M,rit ~ 2.86 X 10-1^ M^'J^M^'^'ptll Mq yr'^ (8) 



We interpret the different masses of the main sequence companion as resulting 
firstly from varying initial masses, and secondly from the star being observed 
at different times in its evolution. 
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6 Population Synthesis 



We calculated the total number of expected X-ray transients evolving through 
case C mass transfer in the same way as Brown, Lee & Bethe (1999). From 
the discussions in Sec. 2, for a ZAMS star of mass 20Mq, we see that only 
the interval of 235i?0 between 137Oi?0 and 1605/?© of a{t = 0) is available 
for case C mass transfer. It is this small fractional interval that allows us to 
obtain a narrow interval of values for Xoce- This is a logarithmic interval of 
only ln(1605/1370) = 0.16 compared with our total logarithmic interval of 
ln(4 X 10^ km/4 x 10^ km) = 7 (Brown, Lee, & Bethe 1999) so that we have a 
fraction of only 0.023 (as compared with 0.11 in Brown, Lee, & Bethe 1999). 

Since the possible ZAMS range for case C mass transfer is very narrow near 
2OM0, we somewhat arbitrarily choose a ZAMS mass interval of I9M0 < 
M < 22Mq. This means (Bethe & Brown 1998) that the fraction of binaries 
with primaries in this range (assuming an IMF, P{m) oc M~^'^) is: (1.9)"^^^ — 
(2.2)~^/^ = 0.08 which is smaller by a factor of 2.5 from that for the interval 
of 20 - 35M0 chosen by Brown, Lee & Bethe (1999). 

From the discussion in Sec. 3 and 4, the upper limit of the mass of the donor 
star is sensitive to the efficiency parameter Aoce- For Aoce = 0.4 this upper 
limit is 1.6 Mq. Assuming a fiat q distribution, we get Ag 1/20. From Sec. 5, 
we take the average lifetime of X-ray transients as 2 x 10^ yrs, a factor two 
higher than the assumed lifetime in Brown, Lee & Bethe (1999). 

By taking the supernova rate as 2 x 10~^ yr^^ per Galaxy, we have the number 
of X-ray transients in our Galaxy as 



2 X 10-2 yr"^ ^ ^ ^ ^ ^ ^'^^ ^ "^'^^^ x 2 x 10^ yr ^ 1840, (9) 

where a binarity of a 1/2 is considered as in Brown, Lee & Bethe (1999). The 
total numbers in the Galaxy of such systems is estimated to be between a few 
hundred and a few thousand (Ergma & Fedorova 1998). Q Our estimate is thus 
consistent with theirs. 

Our birth rate for the black hole binaries obtained from Eq. (9) is 10"^ yr~^, 
roughly half the Portegies Zwart et al. rate for binaries with a neutron star. 



^ Of the Wijers (1996) lower limit of 3000 transient black hole sources, 6 out of 9 
sources had a main sequence companion with a short period, so he would have had 
2000 of the latter. 
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7 Black Hole X-ray Bineiries with (sub) Giant Donors 



Transient black hole binaries with (sub)giant donors might be expected to 
have followed the same scenario, but with initially higher companion masses 
and larger separations after the common envelope, so that they start mass 
transfer only when the donor has evolved off the main sequence. However, 
the large space velocity of Nova Scorpii 1994, best explained with a large 
amount of mass loss (5-IOM0) during the explosion in which the black hole 
was formed (Nelemans et al. 1999), and the high mass of the black hole in 
V404 Cyg (~ IOMq) suggest that the hehum cores of stars around 20 Mq 
may not be massive enough to explain these systems. Hence, the progenitor 
masses of the black hole may have been larger. 

Brown & Lee (2001) require a He core mass of ~ HMq corresponding to 
ZAMS ~ 35M0 for the black hole progenitor of Nova Scorpii and a somewhat 
higher He core mass for the black hole progenitor of V4641 Sgr. The wind 
losses employed by Schaller et al. (1992) must be substantially reduced if 
these SXTs are to be evolved in Case C mass transfer. This would take us 
back to the interval of ZAMS masses 20 — 35Mq suggested by Brown, Lee, & 
Bethe (1999), possibly even up to ~ 4OM0. 



8 Conclusion 



We evolve the low-mass black-hole X-ray binaries, which are identified as the 
observed soft X-ray transients, and show that these systems are formed via 
case C mass transfer following helium core burning phase of the black hole 
progenitor. Although this scenario requires a set of fine tuned conditions, we 
are able to produce a current Galactic population of ~2000 objects, in agree- 
ment with estimates based on observations. Combining the narrow interval in 
initial separations leading to case C mass transfer with the allowed narrow 
range of separations after the common envelope evolution, we put constraints 
on the common envelope efficiency as Actce ~ 0.2 — 0.5. 

Since our analysis requires fine tuned parameter space, more detailed calcu- 
lations of the stellar evolution of stars with ZAMS masses around 2OM0 are 
required. More uncertain is the common envelope efficiency, which is essential 
for the formation of final short orbital period of the binaries. 
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